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ABSTRACT 


Submitted to the Department of Naval Architecture and Marine 
Engineering on May 19, 1967 in partial fulfillment of the require- 
ments for the Master of Seience Desree in Electrical Engineering 
and the Professional Degree, Naval Engineer. 


Electrohydrodynamic power generation has become a field of 
intense interest during the past few years. The utilization of 
electric fields for power generation does permit the use of a lignht~ 
er generator structure which is relatively free of moving parts. To 
date, most of the efforts in this field have been devoted to DC 
power generation. AC power generation is also possible if traveling 
wave synchronous interactions, similar to those ermloyed in the 
conventional synchronous generator, are used. 


An examination of the field equations of the EXD cenerator 
demonstratesthat the exmressions for the develoved power closely 
parallel those of the conventional synchronous machine. Considera- 
tion of the long wavelength limit shows that the power expressions 
can be manipulated in such a fashion as to infer an equivalent 
circuit model of the EHD generator which is the exact dual of the 
cylindrical rotor synchronous generetor equivalent circuit. 


A prototype EHD generator was built to demonstrate the quali- 
tative aspects of the device and to compare the performance of the 
generator with that predicted by the equivalent circuit model. 








age 


Reasonable agreement with the model nredictions was achieved. 
The prototype demonstrated many of the associated problems which 
must be overcome before the device can be optimzed. Careful design 
considerations can eliminate or reduce these prodlems appreciably. 
The AC traveling wave synchronous generator shows promise of being 
able to comrmete with the DC SEHD generator in a relatively short 
period of time. 


Thesis Supervisor: James R. Melcher 
Title: Associate Professor of Electrical Engineering 
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Developed Electromagnetic Power 
Angular frequency (radians per second) 
Rotortii’ wave 
StatorliF wave 
Phase angles 
tream velocity 
Frequency (cycles per second) 
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Wave numoer 
Electric field potential 
Surface charge 
Pereilonsd Electrohydrodynemi.e power 
Permittivity of free space 
linear dimensions 
Power delivered to a resistive load 
Conductance (mhos) 
Resistance (ohms) 
Capacitance (farads) 
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Current source current level 


Average current from current sources 
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Toad current 
Micro (107°) 
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Complex number 
Cormplex conjugate of N 
Capacitive reactance 
Operator "del" 
Electric field 
Stress tensor 
Effieseney 

ateral area 


Particle velocity 
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"Real vart of" 


Time average 





INTRODUCTION 


little attention has been focused in the past on the utilization of 
electric fields for power generation applications. It has long been 
realized that the power densities which can be generated by means of 
electric fields do not compare favorably with those which can be achieved 
by utilizing magnetic field interactions. Large magnetic field densities 
can be achieved through the use of ferromagnetic materials but this 
implies that the generator must inherently possess a high density factor. 
Weight was never a major problem in the past. However, the space ase 
has now imposed the possibility of a major weight problem, and the means 
by which power may be generated through electrohydrodynamic coupling has 
become a significant area of investigation. The major limitations on 
Such devices is the Bilectrical breakdown strength of the environment in 
which they are functioning, and the amount of charge which may be 


deposited on the fluid. 
Magnetic AC Synchronous Interaction 


The majority of standard AC power generating is accomolished through 
ferromagnetic synchronous 3 phase generators. These machines utilize 
the interaction of two sinusoidal m2gnetic fields rotating in Space at a 
synchronous speed. One of the fields, the rotor field, is obtained by 
applying a DC voltage to the rotor of the machine through slip rings. 
The Bemeoiaa) distribution is achieved by the winding connection on the 


rotor. When the rotor is turned by means of a prime mover, the magnetic 





field becomes a sinusoidal rotating (traveling) wave in space. The 
rotating stator magnetic field is also achieved by winding considerations. 
The windinzs on the stator are distributed about the inner periphery and 
are spacially phased at an angle of 360 electrical degrees apart where n 
represents the total phases in the ee winding. As the rotor sweeps 
past the windings, a sinusoidal standing wave voltage is induced in each 
of the phases. This voltage gives rise to a sinusoidal magnetic field 
in each phase, whose amplitude is dependent upon the load demands placed 
on the generator. The resultant wave of the n-phased standing megnetic 
fields is a rotating sinusoidal magnetic field. Thus, the rotor and the 
stator fields both rotate in space at the same synchronous Speed although 
they are spacially separated by an angle ©, (This angle will change 
whenever the magnitude of either the eee field or stator field is 
changed.) 

The expression for the electromechanical power developed by such a 


device can be simplified to the following form: 


Py, = Ky, w FpFS sin @ (1) 


where, w is the synchronous angular frequency 
F,, is the magnitude of the rotor MMF wave. 
ee is the magnitude of the stator *iP wave. 
Kn is a constant reflecting dimensional factors of the 
generator and the parmeability of the sis ZADe 
© is the spacial phase angle by which the rotor field 


and stator field are separated. 





EHD Synchronous Interactions 


Turning our attention to EHD synchronous machines, it can easily 
be demonstrated that marked similarities do exist between the behavioral 
characteristics of the EHD machine and the magnetic AC synchronous 
machinee Consider a linear synchronous machine of the basic form shown 


in Figure l. 


ON § = Re Tes(wt-kx) 
V = Ved6 


Figure 1: A linear Synchronous AC Mechine. 


A very thin stream of charged particles, whose width extends into 
the paper, moving with a velocity U is introduced into a chamber. The 
stream of particles may be described as Of = Re 0 eSut-kx) , where k 
4s the wave number and is equal to 2a » 2 is the wavelength. % is 


the magnitude of the maximum surface chargee The relationship between 


k, Ww, and U is as follows: : 


U=fr= wr=Qy | | (2) 


W, 
out k 


The walls of the chamber have a sinusoidal traveling voltage wave 





impressed on them which has the form V = Re 7 ¢jlut-kx), The e aetes 
wave is moving down the channel with the sams velocity as the stream of 
charged particles. The time dependence of both sinusoidal distributions 
is the same, thus the voltege wave and the stream of charged particles 
are locked in Seni si as they move down the channel. The complex 
voltage V = Vo , 56 where 6 is the spatial phase angle between the 
particle stream and the impressed traveling voltage. Traveling voltage 
waves may be approximated by electrodes, suitably spaced, which have 
n-phased and properly sequenced voltages impressed on them. The SHD 
machine channel is perfectly See with respect to the particle 
stream. 

The solution to the field problem suggested by this configuration 
is developed in Avopendix AI. The final ete shows that the electro- 
hydrodynamic power developed Bee GP being converted to electrical 


energy is: 


Pe = K, w Vo 


O5 sin 6 
This equation may be-modified to the following form by simply 


changing the constant term of the equation: 


Pa = K, w Yo %_ sin & (3) 
d 2€5 


where Vo represents the order of magnitude of the maxamun 
field intensity which can be anticipated in the chamnel if the voltage 


wave alone were present, Ca) and % is the maximum field intensity 
ZEo 
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just above or below the stream of cherged particles due to the presence 
of the charge alone (2c). 


Thus equation (3) may be rewritten in the following form: 
Py’ = K, wo EYE, sin & 


where w is again the synchronous angular frequency and K, is 
a constant reflecting the dimensional factorsof the generator and the 
permittivity of free space. 
Equations (1) end (4) are identical in format. Thus the electric 
and magnetic field synchronous interactions suggest that a possible 


dualism exists in these mechanisms. 


Generator Behevior in the Presence of a 


Contimous Resistive Loading 


The next step in the analysis of an EHD synchronous generator is 
to consider pure generetor action. This entails the introduction of the 
stream of charged particles, but the traveling voltage wave is induced 
on the electrodes by the charge. This is analagous to inducing 2 
voltage in the stator windings of the standard synchronous machine 
through electromagnetic coupling. The locdins of the generator will 
be accomplished by connecting a continuous resistive load to the genera-~ 
tor channel. The load is connected to a ground plane which is held at 
a zero potential. The load is described by a linear conductance of G 
mhos per meter in the X direction. Figure 2 shows a aes ae 


of the generator in the loaded condition. 





jee | / 
| , Aust ioe) 
b region (a) é. = Rega eJ 
1 Be oe. G'(nhos_per meter) _ _ _ _ _ - 
/ > : - 
d region (b) ‘= Regd eJ(ut=1x) 
f° = Re % 5 i(wt=lee) U 
d region (c) é. = Rape ed (wt-kcx) 
a {- ms (ses nee meter) ; 
b region (d) bo _ Regd ed lwt=kx) 


Figure 2: Top View of a Contirmously Loaded AC 


EHD Synchronous Generator 


The field equations of this model are developed in #ppendix “« II. 
The final exoression for the time average power per unit length 
delivered to the load is: 


G)% sinh ko] 
2 Cenk 





<P’ = 


one a 


G \e 
cosh® k(b+d) amen sinh kb* cosh® kd 








Although the equation appears quite complicated, consideration of the 
long wavelength limit which implies that A>? b or d, simplifies the 


equation to the following form: 


ae Bul .-caae) = 
2 Zee. 
8 [w ae ex + G~ b | 


Appendix A III rearranges this equation in such a manner as to suggest 


an equivalent circuit model for the machine which is shotrm in Figure 3. 


€ (a/b) R=1/¢ 





Ficure 3° The Proposed Equivalent Circuit Model for 
the EHD Synchronous Generator. AIL 


quantities are for a per unit length basis. 


The equivalent circuit applies on a per unit length basis. If. the 
unit lencth selected is that of the electrode width, then the circuit 
represents an equivalent circuit model for each pair of electrodes. 


This circuit verifies the dualistic behavior of the two types of 





’ synchronous machines, since the equivalent circuit of the synchron- 
ous ferromagnetic generator is a So liaccustenne in series with an 
inductive reactance and a load imvedance. 

Maximum power transfer in such a circuit is achieved when the 
resistance is equal to the shunt cavacitive reactance. A graphical 
representation of power delivered to the load as a function of load 
resistance is showm in Figure 4. If the load resistances can be 


achieved, the curve can be the basis for checking the validity of 


‘the model. 


~~ —- 





Figure 4: The variation of power delivered to the 
load as a function of load resistance 
for the proposed circuit model. 
Although the entire development of the electrohydrodynariie 
synchronous interaction princivles is based on a continuun medium, 
the Pe aeteeion of discrete electrodes does permit the continuun to 


be approached exverimentally. 





‘DESCRIPTION AND DESIGN HISTORY 


OF THE GENERATOR 


The generztor was primerily a wind tunnel tyve of structure. The 
main air stream velocity (U) was provided by 2 250 CFM axial flow fan 
with a three phase induction motor drive. The air velocity was 
controlled by a three phase rheostat in a tere, line which governed 
the slip of the induction motor. A diffuser section of plexiglass and 
wood concentrated the air stream before entering ae main generator 
section. An aluminum flow straightener was mounted in the diffuser. 
The two remaining sections, the exciter and the electrode chamber, were 
the primary components of the generator. 


2 


The Exciter 


The exciter! proviced the stream of charged particles for the 
Penettitor. The first attempts to investigate a possible exciter used 
a fine tubing mechanism which released particles under a pressure head 
into the air stream. These particles were subjected of a charging 


potential of 2000 to 4000 peak volts at a 60 cycle per second frecuency. 


lone design of the exciter section is the product of a concurrent thesis. 
For a more elaborate discussion of the exciter see TRAVELING WAVE 
CHARGED PARTICLE GENERATION, L. HEROLD, SM Thesis, Massachusetts 


Institute of Technolory (1967). 
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(Sixty cycles per second were used throughout the experiment.) The 
output signal was very random, primarily due to the fact that the release 
of particles in this manner was itself random. Multiple tubes were then 
used, but much of the same effect was noted. Next, commercial spray 
nozzles were tried. These nozzles released a soray of very fine perticles. 
An annular electrode was mounted near the nozzle throat, and with a 
charging potential of only 50 volts RS, a reletively pure sinusoidal 
output was achieved. Peak signal strength was increased by two orders 
of magnitude due to the increase in charge transporting canability which 
is functionally related to the concentrated field strength and the 
increase in particle density. Random signals were no longer encountered. 
The commercial nozzles suffered from the fact that they had to be mounted 
perpendicular to the air flow and a Gant of the particles simply 
Sprayed on the opposite sides of the electrode chamber. A large 
percentage of the charged perticles were lost tefore they effectively 
interacted with any of the electrodes. 

The commercial soray nozzles, hovever, determined the subdsecuent 
design patterns which followed. It was evident that small particles 
in copious quantities were necessary, and that reduced charging 
voltages could accomplish the desired effect easily if proper field 
strengths were used. A series of spray nozzles was then designed which 
introduced fine particles Beets inte the air stream from thin ele- 
ments mounted vertically in the air flow. The final Reser used an 


ejector principle to introduce the particles into the air stream. 





Perpendicular water headers were mounted on the exciter walls. Thin 
tubular connections led from the header to the verticle elements. 
Compressed air was introduced at the top and bottom of a circuler 
section drilled through the center of the vertical element. The air 
was released through small nozzle outlets mounted in the do:mstrean 
face of the vertical elements. The exhausting air crested a vacuum 
in the header and the water was drawn up into the header from a water 
resevoir and forced out through the nozzle openings. The electrode 
for charging the particles was painted on the dowmstream face of the 
vertical element with a silver conducting paint. The upstream face 
was ellipsoidal in shape to reduce air turbulence, while the down- 
stream face was a blunt flat surface. The water connection was used 
for a ground point. A charging potentiel of 200 volts RMS produced 
an increase in peak signal strength of two more orders of magnitude. 
A net overall improvement of five orders of magnitude in the ratio of 
peak output current level to charging potential was achieved in reach- 
ing the final design. Figure 5 shows the exciter section which was 
used in making the load measurements on the generator. Another 


exciter section releasing particles under air pressure is shorm in 


Figure 6- 





Fisure 5: 





The exciter section used for measuring 
the generator performance. Note the 

vertical headers and the connections to. 
the vertical elements. The nozzles are 


mounted in ths face. 








Figure 6: An exciter section delivering a stream of 


water particles. 





The Electrode Chamber 


The exciter section married directly to the electrode chamber. 
The electrodes were painted on the interior of the channel with a 
conducting paint. There were eighteen pairs of electrodes in all. 
The electrodes had similar dimensions, one inch in width and four 
inches in height. Connections led from the electrodes to the exterior 
load resistors. In addition, one megohm resistors are connected in 
series with the load resistances for current measurement purposes. 
Resistance levels in EHD applications are so high that one megohn 
approaches short circuit conditions. 

The electrodes, in addition to being one inch wide, were spaced 
one inch apart. These distances were selected so that the inter- 
action of the particle stream and the traveling wave could be studied 
using impressed high voltages from a six phase transformer. These 
distances would allow a peak potential difference of 30KV bdetween 
electrodes before breakdown would occur. However, this latter group 
of experiments was abandoned vdecause it would add little additional 
insight, and was relatively difficult in comparison to the pure 
generator studies. 

The channel was constructed of 3/8 inch plexiglass and was 36 
inches in length, 2-1/4 inches in width, and had an overall height 
of 6-1/4 inches. These dimensions were the inner dimensions of the 
channel. The contemplated wavelength was to be one foot with six 


electrodes per wavelength. This wavelength would require a lominar 
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flow velocity of 60 feet per second. Although the flow velocity was 
Bemevsble with the fan used, conditions were turdulent rather than 
Jaminar, and the use of compressed air in the exciter did not improve 
conditions either. The compressed air alone yielded 2 flow velocity 
greater than 60 feet per second. The actual wavelengths achieved 
were much greater than one foot, but the ratio of wavelength to the 
width dimension permitted the use of the long wavelength limit in the 
solution to the field equations. Figure 7 shows the entire generator 


| aSSETIDLY « 





Figure 7: The entire generator assemoly. Note the electrodes 


and load resistor connections. 
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RESULTS 


The design of the basic generator was centered around practical 
considerations and the availability of materials. This, however, 
introduced severe limitations. The flow velocity of the fan was in 
a region of marked turbulence. The interaction of the nozzle air 
Supply with the fan was such that control of the resultant velocity 
of the air stream was limited to a semll range. Much higher veloci-~ 
ties are necessary to contain the spread of the particle flow, and 
these higher velocities in turn will more closely approach laminar 
flow conditions within the channel. 

Wetting of the channel proved to be a severe problem also 
because of its effect on the loss of effective charge. Figures 8 
and 9 show the rapid deterioration of the signal in the channel. The 
combination of turbulence and wetting yielded an erratic current level 
pattern but the pattern basically repeated itself under all conditions 
of load. Uniform current levels were finally obtained near the end 
of the channel where turbulence effects can be considered mollified 
and the wetting was severe but uniform throughout. However, the out- 
put current in this region was drastically reduced. 

Another effect noted was the fact that the rate of wetting can be 
increased by increasing the charging voltage on the exciter. This was 


to be expected because of the increasing reoulsive forces between the 
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particles. Once wetting had started, it continued until it reached 
an equiliorium saturation point and little could be done to halt its 
progress with the present design. 

In order to test the validity of the circuit model approach to 
the entire generator, some adjustment in the theory was necessary 
because of the variation in the output current levels. To this end 
the entire generator was modeled as a device cont2ining three bands 
of current generators with six generators per band. The current 
generators in each band possess a current I which is an average of 
the current levels in each of the three bands. Figures 8 and 9 also 
show the band level selected. Passive paramstersfor the equivalent 
circuit were chosen from typical electrode loading behavioral 


¢ 


characteristics for e2ch band level. 


Equivalent Circuit Modeling of the Generator 


The behaviior of electrodes under conditions of increased resis- 
tive loading are shown in figures 10 and'1ll. Tne electrodes snow 
(2, 9, and 15) were typical of the behavior in each band. The 
reactance was computed at two charging voltages and averaged. Only 
a minor discrepancy was noted in each case, and this would be 
consistent with normal data error. 

The manner in which this was done can be ascertained by a quick 


look at the equivalent circuit. 
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Figure 10: ELSCTRODE BEHAVIOR UNDER CONDITIONS 
OF INCREASING LOAD RESISTANCE 
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OF INCREASING LOAD RESISTANCE 





~22— 





tt 
ji 
(eo) 
= 
i 
“| 
ae] 
? 
} 


Xe 





The current at a load resistance of one megohn effectively renra- 
sents a short circuit current which would be the current of the current 
generator or i,.- The current at a load resistance of 100 megonms was 
selected for ip. This was obtained from figures 10 and ll. 

In order to correct for current levels an averase current odased 


on a square root approximation was used for each band, that is 





ip = nba)* 
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Six pair of electrodes: exist in each current band, hence n = 6. 


The results of the calculations are tabulated below: 


—s 


Blectrode Xe (Ma) C(t felectrode) 16 (12 ) Lo (pa) 


(110%) (2007 ) 
2 157 16.9 0557 12Gu7 
9 130 Zoe 0 376 2670 
16 124, 21.5 en: 6398 


In order to test the validity of modeling the generator in this 
fasnion, the power output of the model was cormuted for several values 
of load resistance and compared to the experimental values achieved. 
Figure 12 shows the results of this comparison. A reasonable agree- 
ment is achieved so that some validity appears to exist for the 


derived circuit model. 
Wavelength Measurements 


One ratner obvious fact not mentioned as yet is the phase delay 
which occurs in the induced voltage on the electrodes. This is 
related to the frequency, particle velocity and the distance of the 
electrode pair from a given reference electrode. let electrode pair 
A be the reference electrode, and electrode pair B be another elect- 
rode located a distance L dotmstream from A. The voltages induced 


in A and B would then be: 
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Vey Re Ne e Just 
a Re Vp e just~kL) 


Many important relationships can be deduced from the phase shift 


angle 8 since: 


8 =kL = 2nL = 2nfL = wh (5) 


Sslomts- aoe 


A Vv 





Therefore, the wavelength and the particle velocity may be measured 
directly from phase shift considerations. Since the pressure and 
velocity in the channel were continually changing, the values 
inferred by these measurements were an average value. The accuracy 
of phase snift measurements from an oscilloscope is somewhat limited, 
but simificant changes can yield aeneieGus inSwente 

Wavelength measurements were taken under two loading conditions 
at the entrance and exit ends of the electrode chamber. The results 


are tabulated below: 


Lent (ft) Vp( ft/sec )ent 
R, = O- 1.403 
Ry = Ven e939 80.0 
leat(ft) v,(ft/sec)exit 
Rg = 0 1.29 7709 
Re = WY 1.28 HBAS 
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In order to check the relative magnitude of these quantities, 
manometer readinzs were taken at the beginning and end of the 
channel to measure the flow velocity. These readings indicated an 
entrance velocity of 77 feet per second and an exit velocity of 67 
feet per second. Complete agreement botueen the two measurements 
is out of the question since the manometer is measuring the auetieey 
at a discrete point in a turbulent air flow, while the oscilloscope 
measurements are based on an averaging process. The agreement 
between the velocity figures is satisfactory, however. 

The changes in wavelength do demonstrate the incre2sing inter- 
action of the particles and the induced voltages as the load resist- 
ance is increased. ‘whether or not this is directly related to the 
Slip indicated by the negative stress tensor of Appendix A I is 
questionable at this point, but the ecules do demonstrate that tne 
change in wavelength is greater in the entrance portion of the 


channel where the effective charze 1s proportionately greater. 
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CONCLUSIONS AND KECOLSENDATIONS 


The basic preriise of the thesis was to demonstrate the qualitative 
aspects of an EHD traveling-wave synchronous generator. This has been 
acconplished. The proposed circuit model does seem to predict the 
generetor behavior within tolerable limits. In addition, the phase 
shift measurements do demonstrate the interaction of the charged parti- 
eles with the induced voltages. 

The generator has also clearly shown the conflicting problems which 
are inherent in the design of such a system. Several of these have been 
ignored in the discussion so far, but are mentioned in the following 


list of salient considerations for perfecting the design: 


1. The flow velocity used is in’a region of a turoulent Reynolds 
number. fueh higher velocities will be necessary to counter- 


act the turbulence effect. 


2e The inherent spreading characteristic of a nozzle causes the 
‘electrode channel to wet. This reduces the effective charge 
present in the channel. Higher flow velocities should help 
to reduce this effect, but nozzle devices are necessary for 


providing the copious quantities of charged particles needed, 


3. Data from the exciter indicates that the charge deposited on 
the water droplets is directly proportional to tho charging 
potential in the charging potential regions tested. <Jncreas~ 


ing the exciter voltage increases the explosion tendshey 








lh, 


5e 


Ve 
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of the particle beam and hence promotes wetting. 


External capacitive effects are far greater than anticipated 
from the predicted values. The capacitance is an order of | 
magnitude greater than the preliminary estimates suggested. 
Jt is this capacitive effect which leads to the early loading 
of the generator. Much of this capacitive effect may be in 


the wiring leads, and can be eliminated easily. 


The generator is capable of operating advantageously with 
electrodes connected in parallel if laminar flow conditions 


can be approached and the capacitive effects reduced. 


The feasibility of operating the generator self excited 
seens very plausible as design improvements are incorporated 


into the basic device. 
The dimensions of the electrodes need to be optimized. 


Particle size and mobility need to be studied in much 


greater detail. 


The problems which mist be overcome in order to perfect this type 


of system ere quite numerous and will entail considerable time and 


genius but certainly many improvernents are capable of being incorpora- 


ted into the basic design. 


The relative merits of this prototype snould te studied in cyapati- 


son to recent developments on a DC EHD ballistic gereretor. iaxinun 
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currents achieved on this device have approached 260 p, a with a charg- 
ing potential of 50 KV. The ratio of outout current to charging 
potential in the DC device is 5.2 X 1077. The AC device produces a 
ratio of peak current to peak charging potential of 5.02 X 1079. 
Paralleling of electrodes in the AC device can improve staging outputs 
far beyond this level. Grids must be employed in order to stage the 
DC outout. 

The prototyps has demonstrated that the future does hold some 
promise for AC EHD generators. The problems in perfecting the design 
are recognizable which is the first step in achieving 2 practical 


encineering system. 
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Derivation of the Expression for Developed ZiD Power 


—.... .... oe 
i (a) 
ry : | Or = Re Ooe5 (uxt=lox) oa 
d (b) 


NB = Re YeS(ist-kac) 
V = V,eJ6 


Figure A-l-1. Top view of an AC EHD synchronous generator 


In regions (a) and (b) Maxwell's equations yield the following: 


‘Ox EB =0. , (Gee) 
Therefore, E=-U¢ and (Aiiee) 
Tet = 0 : (4-1-3) 


Equations (A-1-2) and (A-1-3) imply that laPlace's equation must hold 


in these regions. That is 
SY eo: (A~1-4) 
Finally, because of the traveling wave dependence $ rust ce of the form 


t ron : 
p= Re $ (ye Ilwt-kex) , The geometry of the configuration in view of 


LaPlace's equation generates the form of the solution, which is: 


7 Kn “~ 
é. = A cosh ky + B sinh ky (A~1-5) 
ra = € cosh ky +D sinh ky (A~1-6) 


Sufficient boundary conditions exist to evaluate all of the 
constants. They are as follows: 


“~~ 
(1) At the electrode boundary the potential is V. 








Be lan 

g(a) =¥ b(-d) = 7 

E cosh kd + B sinh ka = F (A-1-7) 
€ cosh kd ~ D sinh ka = f (A-1-8) 


(2) The tangential component of the electric field mst be 
contimous at y = 0. 
A=€ (A-1-9) 
(3) The normal component of the electric field at y = 0 mst 


equal e. that is - 0$.,(0) + d$(0) = % 


, by by o 


- kB + «Kh = (A-1-10) 


oo 


The boundary conditions may be summarized as follows, using matrix 


PIO a 1.ON) « 
cosh kd 0 Sinh kd 0 ‘A v 
0 cosh kd @) =sinh kd e C = v 
il =~) jae Oo B 0 
O 0 “1 il DB |°o/ ice 4 


In order to determine the electrohydrodynamic power which is avail- 
able for extraction as electrical power, it is necessary to determine the 
time average retarding stress tensor acting on the stream of charged 


particles. The total stress tensor 


force per unit area on the stream. 


(i ae) (A-1-12) 


1 
Tx 








é 





eam 


aga 2 bd 


Therefore, <i = €6 (E yt 2° Ey ) | 
(y=0) 
, 4h a - Bb 
= €oty (Ey - E,”) : 
“Y (y=0) (A-1-12) 


The time average value of T,, may be exoressed in the following fashion: 


A A¥LA ab : 
<> = 1/20 foe 2 Ee 
ont ae ot as (A~1-13) 


Evaluation of the desired quantities reduces the expression to the 


constants of the field equations. 


7” 7~ “~~ “~ 
Bee jlaid. ES (y=0) = = glkA (Cian) 
Ey = 0g E* G0) = ~3k 
Ey Ga0) =- Dk (A-1-15) 
AE LA ~ 
KT,> = 2 €o RegkA (B - B) (A-1-16) 
Symmetry requires that B= - D. : 
OAK~ 
<Ty> = cohejk“A B (A~1-17) 


Evaluation of the constants from the matrix shows that, 








A=_V + % tanh kd : (A-1-18) 
cosh kd “2k€e9 

B = -% 
mee (A-1-19) 


Substitution of equation (A-1-18) and (A-1-19) into equation (A-1-17) 


yields the desired result. 
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VO. 


x 
<T,>= - Raikv [9 = - k%0% sin 6 (A-1-20) 
2 cosh kd 2 cosh kd 


The time average electrohydrodynarnic developed power capable of 


being extracted from the stream is, 


P, =1<T,>| au (A~1~211) 





where 4 is the total lateral area of the stream. If the 
chamber has a length L, and an effective electrode or stream height a, 
then the power may be expressed in the following forms: 


P,=k%0% = -LaUsin & GE sey 
2 cosh kd 


or 
onVo O L a Usineo (A~1-23) 
ZX cosh kd 
= w Yo 

2 cosh kd , 


where w is the synchronous. fresuency of the system. 


La sin & (A~1~24)) 


Since discrete electrodes sre ceing used, the expression for power 
developed should be corrected to reflect the fact that the interaction 
occurs only when the charged particles are under the influence of the 
electrodese If an area efficiency is defined, then a proper correction 


factor may be applied. 


Ae 
TN, = A leteral Ceiez5) 


Number of pairs of electrodes x electrode area 


N= A jateral 
= Nac = Ne : 
al. J, (A-1-26) 


where N = numer of v-irs of electrodes 


c electrode width 
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Fa = Pana 


Nac Wes sin 6 


cosh kd 


= K, woo sin 6 (ag=27) 
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Derivation of the Output Power Delivered to a 


Continuous resistive Loading 


The system is now depicted as shown in figure A~2~1. The segmented 


ee 


fe b region (a) é. = Roda J fwt-kx) 
XY 
«+f 

a 


G (mhos per meter) 





region (b) $, = Rego @Jiwt-kx) 
Op = Pa% edlut-lex) : 
<< ee 
d region (c) é. = Rage 9 jwt-kx) 
- i ees oe G(mhos per meter) 
b region (d) | é, = Regd 9 dhut~kx) 


Figure A-2-1. Top View of a Continuously Loaded AC SHED 


Synchronous Generator 


electrodes are piserinacee along the X axis and at Y = -2d. The generator 
is again symmetric about the stream of charged particles, so that the 
solutions to the field equations for the upper regions (a) and (b) are 
equally applicable to the lower regions alsoe Tne resistive loading 
connects the electrodes to ground planes which are located at Y=b, 


and Y'= - (b+d). The electrodes have a height "a" which e:ctends into 





the paper. 


The end view of the generator is as shown in figure A-2-2. The 





Figure A-2-2. nd View of the Synchronous Generator. 


resistive load may be assumed to have a conductance of G mnos per meter 
in the X direction. The end view again demonstrates tne symmotry of the 
generator. In view of the symmetry involved, the potentisl may te 


assumed to have a distribution of the form showm in figure 4-2-3. 





/ // 
region (a) 
region (b) 
ae U 
region (c) 
region (d) 


\ \ vo ae 


Ficure A-2-3. Potential Distrioution of a Continuously 


Resistive Load. 





In order to evaluate the entire solution, it is again necessary 


to consider the boundary conditions. These are: 


1. The normal component of the electric field at y = -d must 


be equal to %/fo. That is 








= Oso + ogo = 
by |v =-a) oy [ty=-a) eco mee 
In addition symmetry requires og? = ~04e ’ 
Oy Ky =-d) oy My = -a) 








So that equation (A~-2-1) reduces to 


OS (yey a! (4-2-2) 


Oy “0 


f 


f 
Ze The potentials mist be continuous at the electrodes. 


é2 (0) - fo (oy="0 (4-2-2) 


36 At the ground planes, the potential is defined as the 


point of zero potential. 


J2(b) = 0 (A-2-1) 


4%. Conservation of charge on the electrode surface provides the 
final boundary condition. If we define the charge per unit 
Jength in the X direction 2s qs, then tno normal components 


of the electric field require that 








~hO~ 





ds = ~a€o| O42 -0fe (A-2=5) 
Oy [Cy=0) dy | (y=0) 








Considering the conductance par unit length in the X 
direction across an electrode of heisht a, conservation of 


charge states that 


-Cds = $(0) G. (A-2-6) 
Ot 


Substitution of equation (A-2-5) into (4-2-6) yields the 
final result. 


juwa “ae (0) = Og - = G 220) (A~2-7) 
Oy Oy 


The singularities in the charge 4s seeieet at the boundaries 
require that the solutions be at least piecewise continuous. The 
Peete tine wave dependence of Of requires that the solution be of the 
form ¢ = Re eS eJMut-kx), In addition, La Place's equation holds 


in regions (a) and (b), hence the geomstry suggests the forn of the 


solution. 
g2 = A sinh ky + B cosh ky (A-2-8) 
$° = sinh ky + F cosh ky (A-2-9) 


.The boundary conditions may now be rewritten utilizing the «mom 
form of the solution. Equation (A-2-2) requires that 


2 egk [€ cosh kd - F sinh kd] = - 0; (A-2-19) 








ipa 


Equation (A-2-3) requires that 
B a F = 14 (A=-2~11) 


while equation (4-2-1) states that 
“~~ 2 “a 
A sinh ko + B cosh ko = 0. (A=2~12) 


Finally, the fourth eauation comes from (4-2~7) which 


states that 
juace, k[A-¢) =cB. (A~2=13) 


These results may be surmarized in matrix form as follows: 


0 0 cosh kd = sinh kd A -0,/2¢ 5% 

0 2 0 +1 . |B = 0 
sinh kb cosh ko 0 ' 0 C 0 

1 56 /wae ok =a 0 F 0 


In order to determine the power output per unit lensth, it is 
necessary to consider a slight disression at this point. The vovver 


delivered to the resistive load is Phe 


Py 


= 6 ¢2(0) z (A-2~74) 
On a time average basis this reduces to 


= G $4(0). 42*(0) = G B.B* 7 (A-2-15) 
2 a 





wl! Dee 


Thus the matrix need only be solved for B to determine the output 


power e 


o 


Oo sinn ko 
= 2e_k 
iy Oo (A-2~16) 





fea 
cosh k{btd) - Waeok sinh ko cosh kd 


Thus, the time average power delivered to the load per unit length 


is as follows: 


g |% sinh i? 
<P> = 4 “Tee (A-2-17) 
G 


wee@et i 


cosh® k(btd) + fea sinh® kb cosh? kd 


y 
A brief examination <P)? shows that an optimum value of G exists for 


maximum power transfer, 


watok ( cosh k(btd) ) . (A-2-18) 


Cont 7 sinh ko cosh kd 


Therefore; the maximum power transfer m2y be determined. 


Se sith ke) 
= Got 2¢ 4k 


hy cosn’k(b+d) 


< Paes (A-2-19) 


Although the equations appear to be somewhat complicated at this 


° e 
tes.nin 


point, it is worthy to consider the long wavelensth liamt of them 


+ 3. vem t Wao 


form, ie e« kb and kd << 1. This reduces the complexity of tho Taune 
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matics considerably. Equation (4-2-17), (A-2-18), and (4-2-19) 


reduce to: 


fe 05” NO ae 





Ce = (A-2-20) 
8 [w? ae E+ + G2 ©] 
G =" (Oca cee 
opt abettns (A-2-21) 
[on 
etd = Coca (A~2-22) 
16"e,, 
The total maximum power output of a chamber of length IL will be 
equal to ~ 
<P> a ae Se? (A~2~23) 
Vmx = webb 
26 
However, this expression should also be corrected for the area 
efficiency mentioned in Appendix A Il. 
t 
Se ares = <P) max Na = One Wabe N (A-2-24) 
16 € 
Oo 
8 
<PYmax == OK OO (A-2-25) 
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APPENDIX A III 


DETERMINATION CF AN EQUIVALENT CIRCUIT MODEL 





aS 


Determination of an Equivalent Circuit Model 


Appendix A II developed the expression for the time average power 
delivered to the load per unit length. Maninuletion of this power 
expression will yield correct dimensional quantities for infering an 
eq“uivalent circuit for a per unit length. Equation (A-2-17) will 


serve as a basis for this development. 


06 sinh we? 
» = Zko5 k 


NID 


<P. 








a a 
cosh® k(btd) + ( wae =) sinh® ko cosh® kd 


Considering the long wavelength limit, the equation reduces to the 


following: 


g| gee? 
ape EES , (e522) 
ee (2 b \* 
WES 6 


Multiplying the nunerator by (2 a)* and noting that G = 1/R where 
W.a 
Ris the resistance per unit length in the X direction equation 


(A-}-1) may be manipulated into the desired forn. 


Palme: 
-— 


NO fe 


<P> = 


a 











se lectool 
21 2) Lwegte/o)] - 


<Pi> = ee (A-3-2) 
; Re * aos) 
"SENET LY 


The choice of an appropriate equivalent circuit is governed by 
the principle of duality. The circuit model for the basic AC synch- 
ronous generator utilizing magnetic fields is 2 voltage source in 
series with a synchronous inductive reactance and a load resistance. 
It appears reasonable that an appropriate circuit model for the EHD 
generator utilizing electric fields would be a dual arrangement 
employing 2 current source in parallel with a capacitive synchronous 
susceptance and a load conductance. To this end, consider the 


following circuit configuration and the expression for the pover 





delivered to the load resistance. 





ol) 7-= 





« e e * 

<P> = a Gr iy ) R, Cieee3) 

zd 
$s Sea (A~3-4) 

J 

Pee 

J. 2 2 ‘ R 

< P)> = O22" (te) ae (A-3-5) 


By comparing equation (A-3=5) with (A-2-2), the equivalence of terms 


can be easily recognized. The analogy can be defined as follows: 


Ife 


Oo w a ee 
2 
ég(a/b) = C (A-3-6) 
R = Ry 


A comparison of the units of the terms defined in equation (4-3-6) 


shows that % w a has the dimensions of amperes per unit leneth, while 


is 
and R have the dimension of ohms per unit length. There- 





a 
wW €o(a/b) 


fore, the equivalent circuit defining the time average power delivered 
to a load resistance on a per unit lenzth basis is depicted on the 


next pase. 








Maximum power transfer in this model will occur when the admttances 
are equal which is equivalent to an impedance matching condition. 


mnay is 


Ronee = live (A-3-7) 


This is equivalent to equation (A-2-21) which states that 


Cont =U a) = oe 


In order to obtain some feeling for the order of magnitude involved, 


consider the condition where ab, at a frequency of 60 cps. 


Gopt= fo = 2m x 60 x 1 
36 1 x 10? 
Cont = “3453 x 1077 mhos per meter. 


The standard electrode spacing utilized in the generator was one ingen , 


therefore: 





ahQu 


c = 8.5 x 3072 ares 
opt /electrode eee ewer 
P / a electrode 
a —_ 
electrode 


In addition, the circuit requires that maxinum time average pover 


delivered per unit length be equal to 


KPynay = Sor war 
lone 


which is identical to the maximim power reauirements of equation 


(A=2-22). 
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